Mon. Not. R. Astron. Soc. 000, 3118-3137 (2011) 



Printed 9 September 201 1 



(MN IATeX style file v2.2) 



Morphological properties of z ~ 0.5 absorption-selected galaxies: 
the role of galaxy inclination 

^ ; Glenn G. Kacprzak, 1 * Christopher W. Churchill, 2 Jessica L. Evans, 2 
^ ; Michael T. Murphy, 1 and Charles C. Steidel 3 



Centre for Astrophysics and Supercomputing, Swinburne University of Technology, PO Box 218, Victoria 3122, Australia 
Department of Astronomy, New Mexico State University, Las Cruces, NM 88003 
California Institute of Technology, MS 105-24, Pasadena. CA 91125, USA 



Accepted June 15, 2011 



ABSTRACT 

We have used GIM2D to quantify the morphological properties of 40 intermediate redshift 
Mgn absorption-selected galaxies (0.03 < W r {2196) < 2.9 A), imaged with WFPC-2/HST, 
and compared them to the halo gas properties measured from HIRES/Keck and UVES/VLT 
quasar spectra. We find that as the quasar-galaxy separation, D, increases the Mgn equivalent 
decreases with large scatter, implying that D is not the only physical parameter affecting 
the distribution and quantity of halo gas. Our main result shows that inclination correlates 
with Mgn absorption properties after normalizing out the relationship (and scatter) between 
the absorption properties and D. We find a 4.3<x correlation between W, (2796) and galaxy 
inclination, normalized by impact parameter, i/D. Other measures of absorption optical depth 
also correlate with i/D at greater than 3.2cr significance. Overall, this result suggests that 
Mgn gas has a co-planer geometry, not necessarily disk-like, that is coupled to the galaxy 
inclination. It is plausible that the absorbing gas arises from tidal streams, satellites, filaments, 
etc., which tend to have somewhat co-planer distributions. This result does not support a 
picture in which Mgn absorbers with W,(2796) < 1 A are predominantly produced by star- 
formation driven winds. 

We further find that; (1) Mgn host galaxies have quantitatively similar bulge and disk 
scale length distribution to field galaxies at similar redshifts and have a mean disk and bulge 
scale length of 3.8 kpc and 2.5 kpc, respectively; (2) Galaxy color and luminosity do not 
correlate strongly with absorption properties, implying a lack of a connection between host 
galaxy star formation rates and absorption strength; (3) Parameters such as scale lengths and 
bulge-to-total ratios do not significantly correlate with the absorption parameters, suggesting 
that the absorption is independent of galaxy size or mass. 

Key words: — galaxies: ISM, haloes — quasars: absorption lines. 



1 INTRODUCTION 

Since the first observational evidence associating foreground galax- 
ies with absorption lines detected in the spectra of background 
quasars (Boksenberg & Sargentll 19781) . researchers have strived to 
determine the relationships between absorbing gas found within 
~200 km s _I and a few 100 kpc of their host galaxies. We have 
yet to develop a complete understanding of the many physical 
conditions under which galaxies predominantly produce their own 
metal-enriched "halos" or accrete material in their immediate envi- 
ronment. It is plausible that most galaxies undergo both these pro- 
cesses at some point during their evolution. Assuming that clear ob- 
servational connections between the galaxy properties and the halo 



* gkacprzak@astro.swin.edu.au 



gas properties are present in nature, our collective hope is that we 
can ultimately understand the conditions that drive outflows versus 
the conditions that indicate accretion (such as galaxy star forma- 
tion rates, morphologies, inclinations, etc., and absorption equiva- 
lent widths, kinematics, and chemical and ionization conditions). 

The Mgn AA2796, 2803 doublet is commonly used to trace 
metal-enriched low-ionized gas surrounding galaxies between 
0.1 < z < 2.5. However, understanding the origins of this gas 
is difficult. Mgn absorption traces a wide range o f neutral hydro- 
gen column densities, 10' 6 < N(Hi) < 10 22 cm -2 ^Churchill et all 
l2000al : lRigbv. Charlton. & Churchilll 20021) . which produces a large 
range of detectable Mgn rest-frame equivalent widths, 0.02 < 
W,.(2796) < 10 A. The large range in Hi column densities implies 
that Mgn absorption arises within a large dynamical range of struc- 
tures and environments that contribute to the complex kinematics 
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seen in typical absorption profiles. It may be that lower column 
density systems trace a different population of structures and gas 
producing processes than do higher column density systems. 

The first suggestion that halos observed via Mgii absorp- 
tion may exhibit a dependence on galaxy-quasar separation was 
a ~ 3cr a nti-correlation between W r (2796) and i mpact parame- 
ter, D (e.g. jLanzetta & Bowenlll990l : lsteidellll995l : ICriurcfiill etafl 
l2000bl) . This could be interpreted to imply the universal property of 
a radially decreasing gas density profile in halos. However, given 
the complex velocit y structure of these absorption sys tems, which 
is independent of D, IChurchill. Steidel. & Vogtl Jl996h argued that 
the absorption likely arises from a variety of ongoing dynamical 
events within the galaxy and halo. 

The first hint of a connection between galaxy mass and 
the presence or absence of absorption was deduced from the 
Holmberg-like luminosity scaling bet ween a charac teristic halo ra- 
dius and galaxy K-band luminosity uSteidelll 19950 - Recent stud- 
ies at 0.2 < z < 1 restrict the range of the luminosity power- 
law slope between 0.2 < ft < 0.5, the characteristic Mgn halo ra- 
dius between 90-110 kpc (for an L* B galax y), and the gas absorp- 
tion covering fraction between ~50-90% lITripp & Bowenl 120051 ; 
IChen & Tinker! 120081 : iKacprzak et al] 120081 ; IChen et alj l2010ah . 
However, t he Mgn gas covering f raction may decrease to ~25% 
by z = 0.1 l lBarton & Cooke 20(jl). It remains unclear what phys- 
ical processes replenish the gas reservoirs and populate halos with 
such high coveri ng fractio ns at high impact parameters. 

iBond et alj fcOOlbh analyzed several systems with 
W r (2796) > 1.8 A at z > 1 and concluded that high equivalent 
width absorpt i on sys tems produced by winds is a plausible model. 
iBouche et alj | |2006|) found an anti-correlation between W r (2796) 
and the amplitude of the cross-correlation between luminous 
red galaxies (LRGs) and Mgn absorbers with W,-(2796) > 1 A. 
Since the cross-correlation amplitude is related to the halo mass, 
they inferred that the majority of strong absorption systems 
are not produced by gas that is virialized within galaxy halos. 
They interpreted the W,(2796)-halo-mass anti-correlation as 
evidence that supernova-driven winds become dominant sources 
of stronger Mgn absorbers, i.e., W,-(2796) > 2 A. Further support 
for a wind interpre tation of such strong absorbers was found by 
IZibetti et alj fe007l) by stacking Sloan g, r, and i band images 
toward 700 quasars. Spectral energy distribution fits to their 
stacked images indicate that early-type galaxies are associated 
with weaker absorption than late-type galaxies, which suggests that 
star formation rates (SFRs) may c orrelate with Mgn ab sorption 
strength. Confirming these results, IMenard et alj d2009h stacked 
over 8,500 segments of SDSS quasar spectra where [On] was 
likely to be present at the redshift of the Mgn absorption and 
discovered a highly significant correlation between the W,-(2796) 
and [On] luminosity. They further demonstrated that the distri- 
bution function of W,(2796) naturally reproduces the shape and 
amplitude of the [On] luminosity without any free parameters. 
This implies that Mgn absorption systems trace a significant 
fraction of the global [On] luminosity function. They suggest 
that star-formation driven outflows are the primary mechanism 
responsible for W, (2796) > 1 A absorption systems. 

The outflow scenario is further supported by the 500- 
2000 km s" 1 winds observed in the high equivalent width and 
asymmetric Mgn abso rption profiles seen i n spectra of highly 
star-forming galaxies i Tremonti et al] 20071 ; Weiner et alj 200S 



iMartin & BoucheJ 120091 : iRubin et alj l2009t) . IWeiner et alj 1 120091) 
further demonstrated that W,(2796) and outflow velocity correlate 
with galaxy SFR. 



The evidences discussed above for strong Mgn systems being 
produced by star-formation driven winds primarily uses data that 
have high detection limits such that W,(2796) > 0.5 A and it is 
likely that these conclusions apply to systems with W,-(2796) > 
1-2 A. In a sample of 71 absorbing and non-absorbing galax- 
ies consistin g of predominantly (85%) W r (2796) < 1 A absorp- 
tion systems. IChen et al] l l2010ah do not find a correlation between 
W r (2796) and galaxy color and suggest a lack of a physical con- 
nection between the Mgn halos and recent star formation history of 
the galaxies. Furthermore. IChen et alj ^201 Obi) found that the Mgn 
halo size scales with the stellar mass and the specific star formation 
rate of the host galaxy. They interpret this result as massive galax- 
ies having more extended halos and that Mgn absorption systems 
ari se from infalling clo u ds that fuel star formation. Taken togethe r, 
thelBond etafl fcOOlbh . lBouche et al] j2006h , lzibetti et all d2007h . 
IMenard et alj J2009h . and lChen et alj d2010allbh results may be sug- 
gesting that W,(2796) < 1-1.5 A absorption systems are not wind 
dominated but that larger W,-(2796) systems are wind dominated. 
Studies that incorporate the gas kinematics, galaxy morphologies 
and galaxy orientations relative to the quasar line of sight over the 
full range of W r (2196) would be useful to examine these ideas. 

lKacprzaketalj|20103) compared Mgn absorption and galaxy 
rotation kinematics of 10, z ~ 0.5, W r (2796) < 1.4 A systems 
and found that, in most cases, the absorption was fully to one side 
of the galaxy systemic velocity and usually aligned with one arm 
of the rotation curve. T hese results are con sistent with an earlier 
study of six galaxies bv lSteidel et alj j2002h . However, both stud- 
ies demonstrated that in virtually all cases, a co-rotating halo model 
poorly reproduces the Mgn absorption velocity spread. This im- 
plies that although disk rotation may account for some of the halo 
gas kinematics, other mechanisms must be invoked to account for 
the full velocity spreads. In only 3/17 cases studied in this man- 
ner, all of which had W, (2796) > 1 A, the absorption kinematics 
of th e systems displayed possible signatures of winds or superbub - 
bles jBond et al.l200ld ; lEllison et alj2003l : lKacprzak et alj2010al) . 
Using quasar absorption line m ethods through cosmological simu- 
lations, Kacprz ak et al. I feOlOal) demonstrated that the majority of 
the Mgn absorption arose in an array of structures, such as fila- 
ments and tidal streams, which were all infalling toward the galaxy 
with velocities in the range of the rotation velocity of the simulated 
galaxy. 

IKacprzak et alj d201 lh also directly compared the relative 
Mgn halo gas and host galaxy kinematics for 13 ~L* galaxies at 
z ~ 0.1. They found that these galaxies had low SFRs, low SFRs 
per unit area, and their NaiD (stellar+ISM tracer) and Mgib (stel- 
lar tracer) line ratios implied kinematically quiescent interstellar 
media containing no strong outflowing gas. Given that these host 
galaxies were in isolated e nvironments and given the relative halo- 
gas/galaxy velocity offsets. [Kacprzak et alj d201ll) suggested a sce- 
nario in which the cool halo gas was infalling and providing a gas 
reservoir that could maintain the low levels of star formation within 
the host galaxies. 

The majority of these studies do not incorporate quantita- 
tive morphological and geometric parameters of the host galax- 
ies. Usmg_^round^basedjma^m of Mgn absorbing galax- 
ies, [Steidel]3jcnnson^&Peresog j 19941) determined that although 
most galaxy colors are represented, the absorbing galaxies have an 
ave rage B - K color consistent with that o f a local Sb galaxy (also 
see IZibetti et alj|2007l ; IChen et alj|2010al) . To date, high resolu- 
tion space-based WFPC-2/f/Sr images of Mgn absorption-selected 
galaxies have been used only to qualitatively state that their mor- 
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phologic al types appear similar to those of local spiral and elliptic al 
galaxies dSteidell 19981 IChen et alfeOOlHChen & Lanzettall2003l) . 

In a first effort to quantify morphologies of Mgn absorbers, 
iKacprzak et ail d2007l) used a two-dimensi onal decomposition fit- 
ting program GIM2D dSimard et ai]|2002h to model the physical 
structural parameters of galaxies and compared them to Mgn ab- 
sorption properties. They found a correlation between the galaxy 
morphological asymmetry, normalized by impact parameter, and 
W r (2796). Their correlation increases in significance when strong 
systems [W,-(2796) > 1.4 A] were removed (the greater scatter for 
stronger systems perhaps suggesting a different mechanism giving 
rise to the gas in stronger systems). Furthermore, they reported that 
Mgn absorption-selected galaxies have higher levels of morpholog- 
ical perturbations than found for typical field galaxies. Their cor- 
relation suggests that known processes responsible for populating 
halos with gas, such as satellite mergers and longer range galaxy- 
galaxy interactions, also induce minor perturbations observed in the 
morphologies of the host galaxies. 

The ability to use the WFPC-2/WSr images for more than 
measuring magnitudes or impact parameters, by actually quanti- 
fying the properties of the galaxies, is a useful way of exploring 
the absorber-galaxy connection and may shed light on the differ- 
ences in results found for different equivalent width regimes. In 
an effort to understand the causal connection between absorbing 
ha lo gas and its host gala xy(ies), we expand on our previous work 
of IKacprzak et all d2007l) and apply GIM2D modeling to WFPC- 
2/HST images to quantify addition al morphologica l prope rties of 
known Mgn absorbing galaxies. In IKacprzak et al.l d2007t) we fo- 
cused on morphological asymmetries of absorption-selected galax- 
ies and how they differ from an "ideal" galaxy, which was mod- 
eled by exponential disk and a de Vaucouleurs bulge (n = 4). In 
this paper, we fit a m ore realistic model to 40 galaxies, 34 from 
IKacprzak et alj J2007t) plus an additional 6 galaxies selected from 
the literature, using an exponential disk and a Sersic profile bulge 
(0.2 < n < 4.0) in order to quantify an additional twelve galaxy 
morphological parameters. We also compare the galaxy morpho- 
logical properties to the associated Mgn absorption properties (in- 
cluding kinematics) and explore plausible relationships between 
the absorbing gas and host galaxy properties In § [2] we describe 
our sample and analysis. In § [3] we present the measured prop- 
erties of our sample and identify correlations between galaxy and 
absorption properties. In §[4] we discuss what can be inferred from 
the results. Concluding remarks are offered in §[5] Throughout we 
adopt an Ho = 70 km s~'Mpc~', Qm - 0.3, Aa = 0.7 cosmology. 



2 GALAXY SAMPLE AND DATA ANALYSIS 

We have constructed a sample of 40 galaxies, with spectroscop- 
ically confirmed redshifts 0.3<z<1.0, selected by the presence 
of Mgn absorption in quasar spectra (1^(2796) > 0.03 A). 
We selected 34 galaxies from IKacprzak et alj d2007l ) and 6 ad- 
ditional galaxies found in the literatur J]. The absorption prop- 
erties were measured from HIR ES/Keck dVogt et alTl 19941) and 
UVES/VLT dDekker et alj|2000l) spectra. Galaxy properties were 



1 Three galaxies are not used from the [Kacprzak et aj] 120071) sample: 
Q1127-145G1 and G2 with tabs — 0.312710 were discovered to be in a 
group of at least five galaxies at similar redshifts close to the quasar line-of- 
sight iKacprzak et alj2010bl) . Q2206-199G1 with z abs = 0.751923 is likely 
a star dKacprzaketalj2010al) . 



Table 1. Keck/HIRES and VLT/UVES high resolution quasar observations. 
The table columns are (1) the quasar field, (2) the quasar redshift, (3) the 
instrument used, (4) the observation date(s), and (5) the integration time in 
seconds. 





Z m 


Instru- 




Exposure 


QSO Field 




ment 


Date (UT) 


(seconds) 


Q0002+051 


1.90 


HIRES 


1994 Jul. 05 


2700 


Q01 17+213 


1.49 


HIRES 


1995 Jan. 23 


5400 


Q0229+131 


2.06 


HIRES 


1999 Feb. 08 


3600 


Q0349-146 


0.62 


UVES 


2005 Oct. 11" 


1200 


Q0450-132 


2.25 


HIRES 


1995 Jan. 24 


5400 


Q0454-220 


0.53 


HIRES 


1995 Jan. 22 


5400 


Q0827+243 


0.94 


HIRES 


1998 Feb. 27 


7200 


O0836+1 1 3 

V/UO JV.JT 1 J f 


2.70 


HIRES 


1998 Feb 26 


5400 


Q1038+064 


1.27 


HIRES 


1998 Mar. 01 


7200 


Q1127-145 


1.18 


UVES 




24,900 


Ql 148+387 


1.30 


HIRES 


1995 Jan. 24 


5400 


Q1209+107 


2.19 


UVES 


2003 Mar. 12" 


14,400 


Q 1222+228 


2.05 


HIRES 


1995 Jan. 22 


3600 


Q1241+176 


1.27 


HIRES 


1995 Jan. 22 


2400 


Q 1246-057 


2.24 


HIRES 


1998 Mar. 01 


3600 


Q1317+277 


1.02 


HIRES 


1995 Jan. 23 


3600 


Q1354+195 


0.72 


HIRES 


1995 Jan. 22 


3600 


Q1424-118 


0.81 


UVES 


2005 Jul. 30 


1440 


Q1622+235 


0.93 


UVES 


2003 Jul. 18° 


9800 


Q1623+268 


2.52 


HIRES 


1995 Aug. 20* 


50,360 


Q2128-123 


0.50 


HIRES 


1995 Aug. 20 


2500 


Q2206-199 


2.56 


UVES 


. . /' 


53,503 



" The PID for Q0349- 146 is 076A-0860(A). The Ql 127- 145 quasar spec- 
trum was obtained over multiple nights. The PIDs for this quasar are 67. A- 
0567(A) and 69.A-0371(A). The PID for Q1209+ 107 is 68.A-0170(A). The 
PID for Q1424 - 1 18 is 075A-0841(A). The PID for Q1622 + 235 is 69.A- 
0371(A). The Q2206- 199 quasar spectrum was also obtained over multiple 
nights for PIDs 65.00158(A), 072A-0346(A), and 074.A-0201(A). 
* The Q1623+268 spectrum was obtained over multiple nights; 1995 Aug. 
21, 1995 Sept. 10, 1996 May 25, 2005 May 31, 2005 June 01, 2006 Jul. 02, 
and 2008 Sept. 25. 

measured and modeled from F702W or F814W WFPC-2/HST im- 
ages of the quasar fields. Figure[T]summarizes the absorption selec- 
tion, spectroscopic analysis and morphological fitting technique. 
Note that three absorbers have two galaxies at the same redshifts 
(Q0450-132, Ql 127-145 and Q1623+268). Since it is possible 
that a pair of galaxies can give rise to the absorption and such con- 
ditions can provide further constraints into the galaxy-absorption 
connection, we do not exclude these pairs from our sample. Below 
we describe our data and modeling techniques. 

2.1 Quasar Spectroscopy 

Details of the HIRES/Keck and UVES/VLT quasar observations 
are presented in Table Q] The total integrated exposure times for 
each quasar spectrum ranges from 1200 to 53,303 seconds. The 
HIRES spectrum of Q0836+ 1 13 was reduced using the Mauna Kea 
Echelle Extraction (makefBp ackage and the remaining HIRES 
data were reduced using IRAttj. The UVES data were reduced us- 
ing the ESO pipeline and the custom code UVES Post-Pipeline 

2 http://spider.ipac.caltech.edu/staff/tab/makee 

3 IRAF is written and supported by the IRAF programming group at the 
National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
NOAO is operated by the Association of Universities for Research in As- 
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Table 2. i/ST/WFPC-2 Quasar Field Observations. The table columns are (1) the quasar field, (2) the quasar redshift, (2) the quasar right ascension and 
declination (3), (4) the WFPC-2 filter, (5) the exposure time in seconds, and (6) the proposal identification number and primary investigator of the WFPC-2 
observations. 





QSO - RA 


QSO - DEC 




Exposure 




QSO Field 


(J2000) 


(J2000) 


Filter 


(seconds) 


PID/PI 


QO0O2+051 


00:05:20.216 


+05:24:10.80 


F702W 


4600 


5984/Steidel 


QO 109+200 


01:12:10.187 


+20:20:21.79 


F702W 


1800 


6303/Disney 


Q01 17+213 


01:20:17.200 


+21:33:46.00 


F702W 


2008 


6115/Zuo 


QO 150-202 


01:52:27.300 


-20:01:06.00 


F702W 


5100 


6557/Steidel 


Q0229+131 


02:31:45.894 


+ 13:22:54.72 


F702W 


5000 


6557/Steidel 


Q0235+164 


02:38:38.930 


+ 16:36:59.28 


F702W 


600 


5096/Burbidge 


Q0349-146 


03:51:28.541 


-14:29:08.71 


F702W 


2400 


5949/Lanzetta 


Q0450-132 


04:53:13.556 


-13:05:54.91 


F702W 


2500 


5984/Steidel 


Q0454-220 


04:56:08.900 


-21:59:09.00 


F702W 


1200 


5098/Burbidge 


Q0827+243 


08:30:52.086 


+24:10:59.82 


F702W 


4600 


5984/Steidel 


Q0836+113 


08:39:33.015 


+ 11:12:03.82 


F702W 


5000 


6557/Steidel 


Q 1019+309 


10:22:30.298 


+30:41:05.12 


F702W 


5100 


6557/Steidel 


Q1038+064 


10:41:17.163 


+06:10:16.92 


F702W 


4600 


5984/Steidel 


Ql 127- 145 


11:30:07.053 


-14:49:27.39 


F814W 


4400 


9173/Bechtold 


Ol 148+387 


1 1 '51-29 399 


+38-?5-52 57 


F702W 


4800 


5984/Steidel 


Q 1209+ 107 


12:11:40.600 


+ 10:30:02.00 


F702W 


3600 


5351/Bergeron 


Q 1222+228 


12:25:27.389 


+22:35:12.72 


F702W 


5000 


5984/Steidel 


Q1241 + 176 


12:44:10.826 


+ 17:21:04.52 


F702W 


5000 


6557/Steidel 


Q 1246-057 


12:49:13.800 


-05:59:18.00 


F702W 


4600 


5984/Steidel 


Q1317+277 


13:19:56.316 


+27:28:08.60 


F702W 


4700 


5984/Steidel 


Q1332+552 


13:34:11.700 


+55:01:25.00 


F702W 


2800 


6557/Steidel 


Q 1354+ 195 


13:57:04.437 


+ 19:19:07.37 


F702W 


2400 


5949/Lanzetta 


Q1424-118 


14:27:38.100 


-12:03:50.00 


F702W 


2100 


6619/Lanzetta 


Q1511+103 


15:13:29.319 


+ 10:11:05.53 


F702W 


5000 


6557/Steidel 


Q 1622+235 


16:24:39.090 


+23:45:12.20 


F702W 


24,000 


5304/Steidel 


Q1623+268 


16:25:48.793 


+26:46:58.76 


F702W 


4600 


5984/Steidel 


Q2128-123 


21:31:35.262 


-12:07:04.80 


F702W 


1800 


5143/Mecchetto 


Q2206-199 


22:08:52.000 


-19:43:59.00 


F702W 


5000 


6557/Steidel 



Echelle Reduction (uves pople All the quasar spectra are vac- 
uum and heliocentric velocity corrected. 

The quasar spectra are objectively searched for Mgn doublet 
candidates using a detection significance level of 5 cr for the ,12796 
line, and 3cr for the /12803 line. The spectra have detection lim- 
its of W, (2796) ~ 0.02 A (5<x). Detections an d sig nificance levels 
follow the formalism of Schneid er et al.l i 19931) and lChurchill et all 
dl999h - We define a single Mgn system as absorption occurring 
within < 800 km s , but this definition had no affect on our results. 
A si ngle absorption system m ay have several kinematic subsystems 
(see lChurchill & Vogtll200ll) - absorption regions separated by re- 
gions of no detected absorption. Each subsystem is defined in the 
wavelength regions between where the per pixel equivalent widths 
become consistent with continuum at the lcr level. 

Analysis of the Mgn absorption profiles was performed us- 
ing our own graphic-based interactive software for local contin- 
uum fitting, objective feature identification and measuring absorp- 
tion properties. We compute the equivalent widths, optical depth 
weighted mean redshifts (Zabs), apparent optical depth column den- 
sities (N a [atoms cm -2 ]), and flux decrement weighted velocity cen- 
troids, widths, and asymmetries ((V), W vs , andA [km s -1 ]), directly 
from the normalized pix el counts (see lChurchill et al .1 1 999L [2000al ; 
IChurchill & Vogtll200ll . for the precise definitions of these quanti- 



tronomy (AURA), Inc. under cooperative agreement with the National Sci- 
ence Foundation. 

4 http://astronomy.swin.edu.au/~mmurphy/UVES_popler.html 



ties). Quantities for the full absorption systems are measured be- 
tween the pixel of the most blueward extreme and redward extreme 
subsystems while omitting pixels consistent with continuum (those 
outside/between the kinematic subsystems). 

In many of the Mgn systems, corresponding transitions from 
other elements (such as Mgi ,12853, Fen ,12344, ,12374, ,12383, 
,12587, ,12600, Can ,1,13935, 3970, and Mnn ,12577, ,12594, ,12606) 
are detected, which further validates the identification of a Mgn 
system, but are not required to validate a system. In Figure [2] we 
show an example Mgn system that has corresponding transitions 
from other elements. 

In addition to the aforementioned quantities measured directly 
from the flux decrements, we modeled the absorption profiles using 
Voigt profile least-square fitting. Voigt profile decomposition pro- 
vides a means to model each complex absorption profile as multiple 
individual isothermal "clouds". Each cloud is parametrized by a ve- 
locity center, a column density, and Doppler parameter. The overall 
decomposition provides a model of the number of clouds (N c i) and 
their individual Voigt profile parameters. 

To best constrain the model, the profile fits simultaneously in- 
corporate all transitions associated with the Mgn system. This is 
particularly useful in cases where the Mgn absorption is saturated 
and there is a loss of component structure. Thus, we use the unsatu- 
rated Fen and/or Can and/or Mgi transitions to constrain the num- 
ber of clouds and velocity centers. All transitions have the same 
number of clouds and the velocities of the clouds are tied for all 
transitions. For this work, we assume thermally broadened clouds 
(no turbulent component), so the b parameters of each cloud are 
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HST Image 



Model Image 



Residual Image 
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Figure 1. — The full version of this Figure is in the electronic version of the paper, not in the printed version. In the printed version we present a subset of our 
full sample. — (far-left) The HIRES/Keck or UVES/VLT quasar spectra of the Mgn /12796 absorption feature are shown alongside the associated absorbing 
galaxy on the right. The Mgn /I2796 optical depth weight mean absorption redshift is the zeropoint of the velocity scale. The tick marks indicate the number 
of Voigt profile components and the red curve indicates the fit to the data. We do not have HIRES or UVES data available for six galaxies, (left) WFPC-2/HST 
images of galaxies selected by Mgn absorption. The images are 10 times larger than the 1.5cr isophotal area. — (center) The GIM2D models of the galaxies, 
which provide quantified morphological parameters. A scale of one arcsecond is indicated on each image along with the physical scale computed at the Mgn 
absorption redshift. — (right) The residual images from the models, showing quality of the fit and the underlying structure and morphological perturbations of 
the galaxies. The encircled arrow provides the direction to the quasar (galaxy-quasar orientation). The cardinal directions are also shown and the quasar name 
and redshift of Mgn absorption is provided under each set of galaxy WFPC-2, model and residual image. 
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tied for all transitions (yielding the broadening temperature to be 
the same for all ions). The column densities of transitions belong- 
ing to a given ion are tied, but the column densities of different ions 
are freely fit. 

In Figure [2] we show the least-squares Voigt profile model 
(red) of a saturated Mgn doublet with associated transitions from 
other ions. The number of individual Voigt profile components 
(N c i), or "clouds", are indicated with vertical tick marks. 

The fitting philosophy is to obtain the minimum number of 
clouds that provide a stati stically reasonab le model (i.e.,^,, ==1). 
We use the code MINFIT l lChurchill|[l9971) . which takes a user in- 
put model and performs a least-squares fit while minimizing the 
number of clouds. With each iteration, the least statistically sig- 
nificant cloud (based upon fractional error criterion) is removed 
from the model and a new model (less the one cloud) is fit. An 
F-test is performed to determine whether the two models are sta- 
tistically consistent at the 97% confidence level. If they are consis- 
tent, the cloud is removed from the model. Each cloud not meet- 
ing the fractional error criterion is tested in this way until the all 
clouds meet the criterion and the number of clouds is minimized 
via the F-test iterations. We applied the fractional error criterion 
y/(dN/N) z + (db/bf + (dz/z) 2 < 1.5, where dN/N, db/b, and dz/z 
are the column density, b parameter, and redshift fractional errors 
of the cloud. The resulting models are not sensitive to the input 
model as long as the input model contains a sufficient number of 
clouds. For furthe r description of th e Voigt profile decomposit ion of 
Mgn systems seelChurchilll jl997h . lchurchill & Vogtl l l200lh . and 
IChurchill. Vogt. & Charlton] J200ll) . 

We compute a total system Voigt profile column density, N vr 
[atoms cm -2 ], and its uncertainty by summing the individual cloud 
column densities over the total number of clouds. This number is 
comparable to the apparent optical depth column density in unsat- 
urated systems, but provides a slightly more robust total column 
density in saturated systems where the apparent column density 
provides only a lower limit. 



2.2 HST Imaging of Quasar Fields 

Details of the WFPC-2/HST quasar field observations are presented 
in Table[2] The WFPC-2/HST F702W and F814W images were re- 
duced using the WFPC-2 Associations Science Products Pipeline 
(WASPPlI). WASPP data quality verifications include photometric 
and astrometric accuracy and correctly set zero-points. The F702W 
filter provides a bandpass similar to a rest-frame Johnson S-band 
filter for galaxies at z ~ 0.6. The F814W filter provides a bandpass 
similar to the rest-frame S-band for galaxies at z ~ 0.85. Galaxy 
photometry was performed using th e Source Extractor (Sextrac- 
tor) package l lBertin & Arnoutsll 19961) with a detection criterion of 
1.5 a above background with a minimum object area 75 pixels. The 
i»f 7021V and m F guw magnitudes wer e measured using the WFPC- 
2/HST zero points l l Whitmorel 1 9951) , based upon the Vega system. 

Galaxy absolute 6-band magnitudes, M B , were determined 
from the fc-corrected observed m F -jQ2w or m f8lw . To com- 
pute the fc-correctio ns, we adopt e d spe ctral energy distribution 
(SED) templates of iKinnev etal] dl996h . The adopted SED for 
each galaxy was based upon its rest-frame B - K color ob- 
tained from ISteidel. Dickinson. & Perssonl Jl994h . For galaxies 



5 Developed by the Canadian Astronomy Data Centre (CADC) 
and the Space Telescope-European Coordinating Facility (ST-ECF): 
http://archive. stsci. edu/hst/wfpc2/pipeline. html 
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Figure 2. The saturated Mgn doublet, Fen and Mgi are shown for Zabs = 
1.017038 absorption system seen in the quasar Q2206-199. The Voigt pro- 
file fits (red) and the number of individual Voigt profile components (N c i) 
are indicated with vertical tick marks. The Fen and Mgi profiles are used 
to constrain the kinematics of the saturated Mgn profile and to derive Mgn 
column densities. 
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with no color information, we adopted an Sb SED since this 
is consistent the with average color of a Mgn absorbin g galaxy 
dSteidel. Dickinson. & Perssonlll994l : IZibetti et alj|2007l) . 5-band 
luminositie s were computed using the DEEP2 optimal M* B of 
iFaber et all j2007l Table 2) in the redshift bin appropriate for each 
galaxy: M* ranges from -21.07 «z) = 0.3) to -21.54 «z> = 1.1). 



2.3 GIM2D Galaxy Models 

For each galaxy, the morphological parameters were quantified 
by fitting a two-component (bulge+disk) co-spatial parametric 
model to its two-dimensional s urface brightness d istribution us- 
ing GIM2D (Galaxy IMa ge 2D) llSimard et all 20021). GIM2D uses 
the M etropolis algorithm ^Metropolis et al.l 1953l : ISaha & Williams! 
1 19941) . which does not easily fall prey to local minima, to ex- 
plore the complicated topology of this twelve dimensional param- 
eter space. Once a convergence point is satisfied, the algorithm 
Monte-Carlo samples the region and keeps acceptable parameter 
sets, building up a solution probability distribution, which it uses 
to compute the median of the probability distribution of each free 
parameter and its lcr uncertainties. This technique results in an op- 
timal model parameter set and asymmetric errors. 

For our 40 galaxies, we fit the surface brightness of the disk 
component with an exponent ial profile an d we fit the bulge com- 
ponent with a Sersic profile JSersidlT968h where the Sersic index 
may vary between 0. 2 < n < 4.0. These model fits differ from 
lKacprzaket"Z|j2007l) where we focus was on galaxy asymmetries 
and how absorption-selected galaxies differ from an ideal galaxy 
with a de Vaucouleurs bulge and an exponential disk profile. The 
models have a maximum of twelve free parameters: 

(1.) Galaxy total flux 

(2.) Bulge-to-total fraction (B/T) 

(3.) Bulge semi-major axis effective radius (r b ) 

(4.) Bulge ellipticity (e b = 1 - b/a) 

(5.) Bulge position angle (9 b ) 

(6.) Bulge Sersic index (n) 

(7.) Disk scale length (r f /) 

(8.) Disk inclination (0 

(9.) Disk position angle (9^) 

(10.) Model center sub-pixel offsets (dx) 

(11.) Model center sub-pixel offsets (dy) 

(12.) Background residual level. 

A PSF-deconvolved half-light semi-major axis radius r h is also 
computed for each galaxy by integrating the sum of bulge and disk 
surface brightness profiles out to infinity using the best fitting struc- 
tural parameters. The half-light radius may be unreliable for cases 
where there are large differences between the bulge and disk posi- 
tion angles. Additional information regarding the structural param- 
eters derived by GIM2D can be obtained from lSimard etal.U2002h . 
GIM2D outputs the various scale length parameters in units of pix- 
els and we use the plate scale of the appropriate WFPC-2 chip to 
convert to an angular quantity. We then use the plate scale and the 
angular diameter distance from the adopted cosmology to convert 
them to linear quantities. 

GIM2D uses the isophotal area to extract "portrait size" 
galaxy images from original images with an area 10 times larger 
than the \.5a galaxy isophotal area, which is chosen such that 
an accurate background can be computed by GIM2D. During the 
GIM2D modeling process, the models are convolved with the 
WFPC-2 point spread function, which was modeled at the appro- 



priate locations on the WFPC-2 chip using Tiny Tim dKrist & Hookl 
120041) . 

The GIM2D outputs were manually inspected to see if models 
were realistic representations of the observed galaxies. In addition, 
we used previously measured galaxy rotation curves t o further val- 
idate several galaxy models dKacprzak et alj[2010al : ISteidel et alj 
120021) . In 4/40 cases, we found that the model settled on an un- 
reasonable solution (galaxies Q0454G1, Q1148G1, Q1354G1, and 
Q2206G2). We reran GIM2D for these four systems using different 
random seeds, which resulted in more a realistic and stable model. 
A GIM2D fail-rate of 10% may be less important for large samples 
of galaxies, however with our smaller sample it is crucial to ensure 
that we produce representative models. 

We note that galaxy asymmetries, and structures such as tidal 
tails and bars, may affect the galaxy isophotal shape, thereby pro- 
viding an under/over estimate of the true galaxy shape and orienta- 
tion. It is difficult to model these disturbed systems and they are not 
straightforward to interpret. Thus, in cases such as the interacting 
pair Q0450G1 and Q0450G2, the model inclinations are likely to be 
biased with respect to the real physical inclination of the disk. Thus, 
the model results presented here should be interpreted with some 
caution. However, despite the simplicity of the two-component 
model (i.e., real galaxies may have more than two smooth compo- 
nents, such as spiral arms, bars, Hn regions, etc.), careful analysis 
of the output models and parameters pr ovides useful informa tion 
regarding the complexity of galaxies (see lKacprzak et al]|2007l) . 

The impact parameters, D, and their uncertainties are com- 
puted using a combination of galaxy isophotal centroids deter- 
mined by SExtractor and GIM2D model offsets. The dominant un- 
certainty in D is derived from the pixel offset of the galaxy isopho- 
tal center obtained using SExtractor and the isophotal center of the 
galaxy model determined by GIM2D. This offset is typically about 
0.25 pixels. There is also a ~ 0.05 pixel uncertainty in the position 
of the quasar based upon centroiding errors of unresolved sources 
in the images. 

2.3.1 Galaxy Significance Levels and GIM2D Models 

The WFPC-2 image exposure times for our sample range from 
600-24,000 seconds, with a typical time of ~4700 seconds. It 
may be of concern that the longest and shortest exposures could 
yield different measured morphological parameters for similar type 
galaxies. However, the problem is not so simple since there can 
be several galaxies of interest in each quasar fiel d image having 
different apparent magnitudes. ISimarcT et al. (2002) studied sets of 
GIM2D simulations to characterize the systematic biases and ran- 
dom errors in the galaxy structural measurements. They analyzed 
simulated WFPC-2 F814W images containing ~5500 simulated 
galaxies. Sky photon noise and detector read-out noise were in- 
cluded along with background noise brightness fluctuations con- 
tributed from very faint galaxies below the detection threshold. 
ISimard et all d2002l) determined that the structural parameters can 
be recovered from the simulated images, equivalent to typical 
WFPC-2 exposure times of ~2800-4400 s, and that the uncertain- 
ties only become significant when m Fi n W > 23.5. Only two of 
our galaxies have apparent magnitudes above the Simard limit with 
the remaining below m HST = 22.6. Both of those galaxies have 
m HST = 23.6 and are in a WFPC-2 image that has a typical expo- 
sure time of 5000 seconds. 

Here we examine the significance level (5 L) for each galaxy 
in our sample, which is defined by the ratio of the galaxy mea- 
sured flux and the flux uncertainty (see Table [3j- The SL is then 
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Figure 3. — (a) The apparent magnitude as a function of exposure time 
for the 40 Mgn absorption-selected galaxies in our sample. — (b) Galaxy 
significance level (SL) as a function of exposure time. — (c) Galaxy SL 
distribution. — (d) Galaxy 5 L as a function of galaxy apparent magnitude. 
Note that the galaxy S L is only dependent rriHST and not the WFPC-2 ex- 
posure time. 



equivalent to an average signal-to-noise ratio per pixel, with higher 
weighting towards the bright pixels. In Figure[3]:(, we show the dis- 
tribution of galaxy apparent magnitudes as a function of exposure 
time. All exposure times cover a wide range of apparent magni- 
tudes. We see no trends with exposure time. 

In Figure (3j>, we show the distribution of S L as a function of 
exposure time. There is a wide range of SL for each exposure time. 
In Figure[3];', we show the binned distribution of SL. The majority 
of galaxies lie within the main peak around SL ~ 80 with others, 
not necessarily from long exposures, at higher 5 L. In Figure[3]:/, we 
show the S L as a function of galaxy apparent magnitude. Note the 
strong anti-correlation. The anti-correlation exhibits some disper- 
sion which can be attributed to background Poisson noise and also 
some variance in the galaxy surface brightness for a given magni- 
tude. 

The strong anti-correlation between SL and m H sT suggests 
that the galaxy S L primarily depends on its apparent magnitude and 
not on the exposure time. Therefore, degrading individual quasar 
field images to a lower signal-to-noise ratio is not necessary for 
the galaxy magnitude range of interest (m HST < 23.5). We dis- 
cuss additi onal tests regarding S/ N and quasar PSF near absorbing 
galaxies in I Kacprzak et alj J2007l> . 



3 RESULTS 
3.1 The Sample 

We have constructed a sample of 40 Mgn absorption-selected 
galaxies between 0.3<z<1.0. In Figure [T] we show WFPC-2 por- 
traits of the galaxies. Note that the galaxy images are 10 times 
larger than the 1.5cr galaxy isophotal area, giving the galaxies the 
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Figure 4. — (a) The normalized distributions of the bulge effective radii 
derived for 120 galaxies from DEEP survey imaged with the F606W filter 
using a fixed Sersic index of n = 4 (blue), 120 galaxies from DEEP survey 
imaged with the F814W filter also using n = 4 (red), our sample having free 
Sersic index fits (solid black), and our sample fit with n = 4 (dotted black). 
The mean values are also shown and the error bars represent the average 
uncertainty from the model fits of all galaxies. All samples have similar 
bulge effective radii distributions. — (b) same as (a) except plotted for the 
disk scale length distributions. Samples have similar distributions except 
our sample has larger disk scale lengths. Note that when the Sersic index is 
fixed at n = 4, the disk scale length distribution shifts towards smaller sizes 
and becomes more consistent with the DEEP sample. — (c) the redshift 
distributions for the 40 galaxies from our sample. The mean redshift value 
is indicated. The 120 DEEP galaxies have been selected to have the same 
redshift distribution as our sample. 
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Table 3. The observed galaxy properties. The columns are (1) the quasar name and galaxy ID number, (2) the optical depth weighted mean Mgii /12796 
absorption redshift, (3) the Mgii /12796 rest-frame equivalent width, (4) the source of the adopted equivalent width measurements, (5) the adopted galaxy 
redshift, (6) the source(s) of galaxy spectroscopic redshift, (7) the projected galaxy-quasar separation (impact parameter), (8) the WFPC-2/HST F702W or 
F814W galaxy apparent magnitude, (9) galaxy significance level in the WFPC-2 images, defined to be the ratio of the measured galaxy flux and the flux 
uncertainty, (10) the B - K colors, (11) the rest-frame B-band galaxy absolute magnitude, and (12) the galaxy B-band luminosity. 
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appearance of all being of similar size; the angular and physical 
scales of each galaxy image is shown. The encircled arrow next 
to the portrait provides the direction to the quasar relative to the 
galaxy. In Table [3] we list the observed galaxy properties for the 
sample, which were mostly derived from the WFPC-2 images. The 
galaxies range between -18.6 < M B < -23.8 and have a impact 
parameter range of 12 < £) < 115 kpc. 

Also in Figure Q] we show the fitted Mgii A2196 absorption 
profiles. The number of clouds is indicated by the tickmarks above 
the fit. We do not have in hand the high resolution spectra of six 



Mgii systems so we can only use their equivalent widths in our 
analysis. Note the wide variety of Mgii /12796 profile velocity 
widths and structures. The Mgii rest-frame equivalent widths range 
between 0.03 < W r (2796) < 2.9 A. In Table[t]we present the Mgii 
absorption properties measured from the profiles and obtained from 
the Voigt profiles fits. 

Along with the HST images shown in Figure Q] we also in- 
clude the galaxy GIM2D model and residual images. The model 
and residual image is displayed using the same dynamic range as 
the WFPC-2 image. As mentioned above, the smooth GIM2D mod- 
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Table 4. Summary of the Mgn /12796 absorption properties. The table columns are (1) the quasar name and galaxy ID number, (2) the optical depth weighted 
mean Mgn /12796 absorption redshift, (3) the rest-frame equivalent width W, (2796), (4) the doublet ratio (DR), (5) the number of clouds (N c i), (6) the Voigt 
profiles fitted system column density (N v „), (7) the Mgn optical depth, (8) The AOD derived column density (N a ), (9) the mean velocity ((V)), (10) the velocity 
spread, (W vs ), and (9) the velocity asymmetry (A). 
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a In cases where the Mgii /12796 absorption is saturated, we use the Mgii /12803 absorption profile to compute log(/V a ) and t. The Mgii /12796 and Mgii 

/12803 absorption profiles are both saturated in only six cases. 
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Table 5. The Mgii absorbing galaxy morphological properties, which were obtained from GIM2D modeling of the WFPC-2///ST images. The table columns 
are (1) the quasar name and galaxy ID number, (2) the optical depth weighted mean Mgii /12796 absorption redshift (z. a bs), (?) the Mgii /12796 rest-frame 
equivalent width W r (2796), (4) the bulge-to-total fraction (B/T), (5) the bulge Sersic index (n) (6) the bulge semi-major axis effective radius (rt), (7) the bulge 
ellipticity (ey = 1 — b/a), (8) the bulge position angle (8b), (9) the semi-major axis disk scale length (f>/), (10) the disk inclination (;'), (11) the disk position 
angle (8j), (12) the galaxy half light radius (r;,), and (13) the^ 2 per degree of freedom for the fit. 
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els do not include degrees of freedom to fit spiral arms, Hn regions 
etc. seen in some galaxy images. For example, one can clearly see 
these structures in the residual images of Q0229+131G1 (spiral 
arms) and Q0454-220G1 (Hn regions). 

In Table [5] we list the galaxy morphological properties ex- 
tracted from the GIM2D models. The 40 galaxies have a mean disk 
scale length of 3.8 kpc, which is larger than the estim ated value for 
the Milky Way of 2.3 ± 0.6 kpc dHammer et al]l2007t) . The sample 
mean bulge scale length is 2.5 kpc, whi ch is comp arable to esti- 
mates of the Milky Way bulge of 2 kpc dRichll 19981) . The average 
Sersic index for the sample is (n)=1.9. We find 23 galaxies having 
well-modeled disk-like bulges with Sersic indices of n < 2 and 17 
galaxies having steeper bulge profiles with n > 2. 



3.1.1 Mgn Absorbers Compared to DEEP Galaxies 

In order to determine if Mgn absorption-selected galaxies differ 
from field galaxies, we compare the bulge and disk size distribu- 
tion from our sample to galaxies having similar properties obtained 
from the Deep Extragalactic Exploratory Probe (DEEP) Groth Strir 



trom the Deep bxtragalactic exploratory Probe (DbhP) uroth strip 
Galaxy Redshift Su r vey (seelSimard et alj|2002l ; IVogt et al.ll2005l ; 
IWeineretal.ll2005l) . ISimard et al.l j2002h used GIM2D to model 
WFPC-2/ftST F606W and 814W images of 7450 galaxies from 
the DEEP survey. We selected a subset of 429 galaxies that have 
both a redshift and magnitude range similar to those represented 
by our sample (see Tabled- The mean redshift for the 429 DEEP 
galaxies is (z) = 0.68, peaking at z ~ 0.9. In Figure [4J; we show 
the redshift distribution for our sample having a mean redshift of 
(z) = 0.51, peaking at z ~ 0.4. Since the redshift distributions for 
the two samples are quite different, we used our redshift distribu- 
tion as a selection function to randomly select a subset of the 429 
DEEP galaxies to make a fair comparison. We randomly selected 
120 DEEP galaxies that reproduce the redshift distribution of our 
sample seen in Figure[4j;- We produced several DEEP sub-samples 
all resulting in the similar distributions of bulge effective radii and 
disk scale lengths. 

The galaxies from the DEEP sample were fit using GIM2D 
where the bulge Sersic index is limited to n = 4, whereas we al- 
lowed the bulge Sersic index to vary between 0.2 < n < 4.0 for our 
galaxy model fits. In order to make a direct comparison with the 
DEEP galaxies, we refit our 40 galaxies with a fixed Sersic index 
of n = 4. It is also worth noting that the majority of our galaxies 
were observed using the WFPC-2 F702W filter, whereas the DEEP 
survey used two different WFPC-2 filters; a bluer F606W filter and 
a redder F814W filter. The model fits to the DEEP F606W and 
F814W images produce similar size distributions (see Fig|4ji,b). 
Given that our galaxy images were taken primarily using a filter 
with a central wavelength between the central wavelength of the 
two filters used in DEEP, we do not expect that the difference in 
filter band-pass used will produce any shifts/differences in the disk 
and bulge size distributions. 

In Figure |4ji, we show the normalized distribution of bulge 
effective radii averaged over several random sub-samples of 120 
galaxies derived from the DEEP F606W images, the DEEP F814W 
images, our sample of 40 galaxies having free Sersic index fit, and 
our sample fit with n = 4. The mean values along with the mean 
fit errors are shown. For our sample, the mean value is consistent 
with those derived from the DEEP sample, although more galaxies 
are found with smaller bulge effective radii in the latter. The K- 
S probability that the DEEP sample and our sample (fit with n = 
4) are drawn from the same population is P(KS) = 0.27, which 
suggests that the two distributions are similar and only differ at 85% 



c.l. (1.1 ltr). Note that for our sample, when the Sersic index is fixed 
at n = 4, the bulge sizes tend to increase on average. As the bulge 
contribution increases, it is reflected as a smaller size distribution 
of the galaxy disks as seen in Figure|4j5. 

In Figure [4jj, we show the disk scale length distribution. Here 
we find that the averaged sub-samples of 120 DEEP galaxies have 
a lower mean rj than found for our sample. On average, our sample 
has more galaxies with larger ra compared to those from the DEEP 
survey. The K-S probability that the DEEP sample and our sample 
(fit with n = 4) are drawn from the same distribution is P(KS) = 
0.27, which suggests that the two distributions are similar and only 
differ at the 73% c.l. (l.llcr). However, if we compare our sample, 
when using a free Sersic index fit, this increases to 3.7cr and the 
mean value is now inconsistent with those of the DEEP survey. 

We find our sample, when fit with n = 4, similar to the galaxies 
in the DEEP sample: differences arise only when we have a free 
index fit. This difference is only due to the fitting and not a physical 
difference. 



3.2 Galaxy Type and Absorption 

It is difficult to classify the morphological type (i.e., E4, Sb, etc) of 
a galaxy using only the quantified morphological parameters pre- 
sented here. While early-type galaxies are expected to have high 
Sers ic indices of n > 2, o ver 25% of spirals have bulges with n > 2 
(see IWeinzirl et all 12009)) . Thus, the Sersic index is not ideal for 
classifying galaxy morphology. Instead we have used the bulge-to- 
total ratio to separate early-type (B/T > 0.5) and late-type galaxies 
(B/T < 0.5). We expect less th an 8% (±3 galaxies) contamination 
by applying a B/T = 0.5 cutoff l lWeinzirl et all2009l) . 

For our sample, we find 27 late-type galaxies and 13 early- 
type galaxies; these numbers are consistent with the distribu- 
tion of galaxy morphologies found in the fi eld environment 
( Ivan den Ber gh 2004: Ham mer et al 1l2005l.l2009l) . In Figure[5ji, wc 
show the distribution of galaxy type as a function of W,-(2796). The 
mean W,(2196) for the early-type galaxies is 1.0 A and is slightly 
lower for late-type galaxies at 0.7 A; both types are accompanied 
by a large spread in equivalent width. 

In Figure [5p, we show the W, (2796) as a function of impact 
parameter, where red points are early-type galaxies and blue points 
are late- type galaxies. The mean impact parameter for each mor- 
phology type is the same: 57 kpc. There are some interesting trends 
to note. For D < 40 kpc, we find a higher proportion of late-type 
galaxies (12/15) that span a large equivalent width range. Beyond 
40 kpc, we find more of a mixture of early and late-types, with late- 
type galaxies dominating (14/25). In this impact parameter regime, 
we also see that late-type galaxies tend to have weaker absorption 
than early-types. The data suggests that late-type galaxies are the 
dominate galaxy type selected by Mgn ab sorption at all impact pa- 
rameters. We note that lChen et alj ( feOlOal) also find more spirals in 
close proximity to the quasar line of sight and find a factor of two 
more early-type galaxies at higher impact parameters (D > 43 kpc). 
A K-S test yields that the equivalent width distributions of the late- 
type and early-type galaxies in our sample differ by 2.1cr. Given 
the low number statistics of our sample and the larger scatter of 
W, (2796) associated with late-type galaxies, we caution any against 
any definitive conclusions being drawn from Figure|5] 

It could be possible that the physical scale lengths, i.e., r d , r h , 
etc., of a galaxy may be related to its halo gas absorption proper- 
ties. Having quantified the morphological parameters of the galax- 
ies in our sample, we have applied generalized Kendall and Spear- 
man rank correlation tests, which accounts for measured limits in 
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Figure 5. — (a) The distribution of the Mgn /12796 rest-frame equivalent 
width, W, (2796), as a function of galaxy morphology. Galaxies have been 
separated into to classes using the bulge-to-total ratio: late-type (B/T > 0.5) 
and early-type (B/T < 0.5). Note that a large fraction of late-type galaxies 
are associated with low equivalent width systems. — (b) The impact pa- 
rameter as a function of the W,-(2796) is shown for all 40 galaxies in our 
sample. The blue points are late-type galaxies and the red points are early- 
type galaxies. We find a 3.2<x anti-correlation. Note the large scatter in the 
correlation and the lack for high D galaxies with high W, (2796). 



we note that, given the classical methods used to identify host 
galaxies (find absorption then search for host galaxy), our sam- 
ple and previous samples are relatively biased toward galaxies in 
close proximity to the quasar, which may artificiall y lead to the ob- 
served anti-correlation. As u nbiased surveys (e.g.. lTripp & Bowenl 
l2005tlBarton & Cookd2009l) are enlarged, we will be able to obtain 
more robust data to examine the significance of the D and W, (2796) 
anti-correlation. Taken at face value, the observed anti-correlation 
between W r (2796) and D would suggest that either the column 
density and/or velocity spread of the absorbing gas decreases with 
galactocentric distance. 

In Table [6] we present the results of additional correlation 
tests between the absorption properties and D. What is immedi- 
ately clear is that all our measures of the gas absorption strength 
decrease with D, although we note that the doublet ratio increases 
with D as expected. VK,(2803), which is less sensitive to satura- 
tion than is W, (2796), anti-correlates significantly with D. Though 
NJ, N v p and t do not show a significant anti-correlation with D, 
they do follow the trend of decreasing with D. Taken together, the 
anti-correlation trends show a consistent signature that, on average, 
Mgn absorbing gas becomes optically thinner with increasing pro- 
jected galactocentric distance from the host galaxy. 

As stated above, the decreasing equivalent widths indicate that 
either the column density and/or velocity spread of the absorbing 
gas decreases with galactocentric distance. However, the facts that 

(1) W,(2803) shows a slightly higher significance than W,-(2796), 

(2) the doublet ratio increases with D, and (3) the mean velocity, 
W vs , and the velocity asymmetry. A, scatters with D (shows no sig- 
nature of decreasing with D), together suggest that a decrease in 
column density with D, rather than just a decrease in velocity struc- 
ture, is likely the dominant behavior of the halo gas. 



the sample dFeigelson & Nelsodll985» . between the galaxy scale 
lengths, Sersic index, n, and bulge-to-total ratio, B/T and the mea- 
sured absorption properties. We find no evidence for trends with 
these morphological parameters (all significance levels are less than 
2cr). This may suggest that galaxy mass, size, and/or global stellar 
distribution is/are not important factor(s) in dictating the absorption 
properties of gaseous halos. 



3.3 Impact Parameter and Absorption 



For many studies an anti-correlation between W,(2796) and 
D has been noted and employed to understand the geometry, 
mean radial density pro file, and cross-sectional covering fraction 



of Mgn halo gas (e.g., Lanzetta & Bowenl 19901; Steidell 


1995; 


Churchill, Steidel. & VoetT 19961: Mo & Miralda-Escudd 


1996; 


Boucheetal. 200€ 


; Chelouche et al. 120081: IChen & Tinketl 


200S; 


Chen et alj 2010al: 


Steidel et alj|201fj|). For our sample, we 


find a 



3.2cr anti-correlation (see Figure[5p) . Though it is clear that pro- 
jected galactocentric distance is a predominant parameter govern- 
ing absorption strength, there is significant scatter in the relation, 
showing that distance from the host galaxy may not be the only 
physical parameter governing the distribution and observed quan- 
tity of halo gas. For example, at various levels of influence, star 
formation rate, environment, or the orientation of the line of sight 
(galaxy inclination and/or position angle relative to the quasar) 
could also play a role. 

As a caveat with regards to interpreting this anti-correlation, 



3.4 Galaxy Color and Absorption 

Given the results o f IZibetti et all d2007h . who statistically show that 
larger W, (2796) is a ssociated with bluer galaxies, and the results of 
iMenard et all j2009h . who find a highly significant correlation be- 
tween [On] luminosity and W,(2796), we might expect to find that 
bluer galaxies in our sample are associated with larger W, (2796). 
We thus tested for a correlations between B - K galaxy color and 
absorption properties in our sample. 

At the range of significance level ~2-2.2c\ dominated by N c u 
N vp , and (V), we find larger absorption quantities tend to be as- 
sociated with the bluer galaxies. The Spearman Kendall results 
are presented in Table [6] If we were to interpret these trends, we 
would infer that redder (early-type) galaxies are associated with 
low column density absorption with fewer kinematic components 
than bluer (late-type) galaxies. We do not find evidence supporting 
a trend of increasin g W r (2796) with decreasing B - K, as found by 
IZibetti et alj d2007h . 

The result fo r our sample is consistent with the results of 
IChen et ail d2010ah . who examined a sample of 71 systems. They 
find no significant correlation between B - R galaxy color and 
WV(279 6). It remain s difficult to directly compare our work and the 
work of IChen et alj because the works employ d ifferent color s. It 
is also difficult to compare our work to those of IZibetti~et al. due 
to the stacking nature of their data analysis, the difference in the 
W r (2796) distribution of the two samples, and the use of different 
colors in both samples. We further discuss these samples and their 
results in Section|4] 
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Table 6. Selected results of the Kendall and Spearman rank correlation tests between the Mgn absorption properties and the associated galaxy properties. In 
column (1) are the tested properties and (2) has the number of systems tested. In columns (3)-(5) are the Spearman correlation coefficient, r s , the probability, 
P s , that the tested data are consistent with the null hypothesis of no correlation, and the number of standard deviations, N s where the ranks are drawn from a 
normal distribution. In columns (6)-(8) are the Kendall r K , the probability, P K , and N Ktr . 
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1.84xl0~ 5 


4.28 


i/D vs. N a 


33 


0.65 


0.0002 


3.70 


0.93 


3.79xl0" 5 


4.12 


i/D vs. r 


33 


0.61 


0.0006 


3.44 


0.87 


0.0001 


3.86 


i/D vs. Nd 


33 


0.58 


0.0010 


3.30 


0.85 


0.0004 


3.54 


i/D vs. N vp 


33 


0.53 


0.0029 


2.98 


0.78 


0.0015 


3.18 


i/D vs. DR 


39 


-0.41 


0.0115 


2.53 


-0.58 


0.0093 


2.60 


i/D vs. W vs 


33 


0.40 


0.0234 


2.27 


0.56 


0.0227 


2.28 


i/D vs. A 


33 


-0.08 


0.6357 


0.47 


-0.13 


0.6090 


0.51 


i/D vs. (V) 


33 


0.03 


0.8552 


0.18 


-0.01 


0.9629 


0.05 
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3.5 Galaxy Luminosity and Absorption 

A Holmberg-like luminosity scaling of the radius of Mgn ab- 
sorbi ng halo sizes has been discussed a t length in the litera- 
ture dSteidellll995l : Ichen & Tinken 120081 : IKacprzak et all 120081 ; 
Chen et all l2010al). For the sample of galaxies presented here, 
Kacprzak et al, 1 120081) showed that the Mgn absorption halo size 
exhibits a scaling with the galaxy luminosity (albiet not with a clean 
cut-off size). Do the absorption properties also exhibit a scaling 
with luminosity? To investigate this question, we examined if there 
is a correlation between absorption properties and M B (effectively 
the galaxy luminosity). We find no evidence of a correlation with 
M B for any the absorption properties. The Spearman and Kendall 
results are presented in Table [6] This would suggest that galaxy 
luminosity is not a predominant governing factor in determining 
the absorption strength, optical thickness, or velocity spread of the 
Mgn gas in our sample. 



3.6 Galaxy Orientation and Absorption 

lLanzetta & BowerJ d 19921) and ICharlton & Churchilll jl996h 
showed that if there are preferred systematic kinematics and spatial 
distributions of Mgn absorbing gas relative to the host galaxies, 
then the absorption strengths and kinematics would be expected to 
follow a predictable behavior as a function of galaxy orientation 
and impact parameter. By orientation, we refer to the galaxy 
inclination, i, and the position of the quasar line of sight relative to 
the major axis of the galaxy, 8, which ranges from 8 = 90° when 
the quasar aligns with the major galaxy axis to 8 = 0° when the 
quas ar aligns with the galaxy min or axis. 

ICharlton & Churchilll i ll 9961) further demonstrated that Mgn 
halos could be disk-like since, statistically, the Mgn redshift path 
density could be satisfactorily explained by extended disks as well 
as by spherical halos. They predicted that if absorbers are disk- 
like in origin, then (1) the mean W r (2796) for a given inclina- 
tion increases with i, (2) the relative number of stronger absorbers 
peaks when the quasar probes the major axis (8 = 0°), and (3) the 
mean ratio of the absorption kinematic spread to the Tully Fisher 
velocity of the galaxy increases with i. Maki ng a direct compari- 
son of galaxy and Mgn absorption kinematics. lSteidel et alj 120021) 
showed that galaxy halos are partially understood by lagging disk- 
like kinematics . Binni ng Mgn absorption profiles by inclination, 
IKacprzak etliflteOlOiJ) found that the mean optical depth and mean 
velocity spread of the absorbing gas was larger for higher inclina- 
tion galaxies. 

From the above predictions, if flattened halos that are co- 
planer with the galaxy disk dominate, we expect the absorption 
properties to correlate with inclination and anti-correlate with po- 
sition angle (be strongest at i = 90° and 8 = 0°). On the other 
hand, it is also possible that star-formation driven winds could be 
a dominant source of ob served Mgn absorption (especially for the 
stronger absorbers, e.g.. lZibetti et al ,1120071 ; I Weiner et alj|2009l) . In 
contrast to AGN driven winds, star-formation driven winds are spa- 
tially distri buted geometrically perpendicular to the plane of the 
galaxy (see lVeilleux et al.l2005l) . Thus, for winds, we would expect 
the absorption properties and velocity widths to be anti-correlated 
with inclination and correlate with position angle (be strongest at 
i = 0° and 6 = 90°). In reality, both scenarios, possibly further com- 
plicated by additional scenarios such as mergers, IGM filament ac- 
cretion, etc., contribute to the presence of Mgn absorption. Are the 
predicted observational signatures of one of these scenarios man- 



ifest in the data, which would suggest predominance of the sce- 
nario? 

GIM2D yields two position angles - one for the bulge and 
one for the disk. We employ the bulge-to-total ratio for adopt- 
ing the best representative galaxy position angle, such that when 
B/T > 0.5 we use bulge, and when B/T < 0.5 we use the 
disk, 9 d . For our sample, we find that the distribution of inclina- 
tions and position angles of galaxies selected by Mgn absorption 
are not inconsistent with having been drawn from a population of 
field galaxies (for which the distribution of inclinations is oc sin 2 i 
and position angles are random). A K-S test for the inclination dis- 
tribution yields a K-S statistic of 0.16 and P(KS) = 0.27 and for the 
position angle distribution yields 0.16 with P(KS) = 0.28. At face 
value, this would suggest that galaxies chosen by known absorp- 
tion are distributed on the sky no differently than galaxies selected 
at random. This might further suggest that halos have a uniform 
spatial distribution around their host galaxies. However, there is a 
wide range of absorption properties (which are governed by the col- 
umn densities and kinematics, i.e., velocity widths) and the overall 
distribution of orientations relative to the line of sight cannot, by 
themselves, discriminate whether absorption properties vary in a 
systematic fashion with orientation. 

We find no correlation between 8 and any of absorption prop- 
erties (the highest significance is 1.4cr); the absorption properties 
are consistent with a random distribution as a function of position 
angle (these results are independent of our B/T selection criteria). 

The results of correlation tests between the absorption proper- 
ties and inclination are similar to the tests with impact parameter; 
we find several trends toward a correlation above the 2. Oct level. 
The Spearman and Kendall quantities are presented in Table|6] The 
strongest trends are N c u N a , and N vp all showing that these quan- 
tities tend to increase with increasing inclination. Taken together, 
these trends could suggest a potential correlation between disk in- 
clination and the column density, number of clouds, and W, (2796). 

Overall, these results are suggestive (not definitive) that the 
observational signature from the predictions of a co-planer geomet- 
ric model dominate over the other plausible scenarios for the origin 
and loc ation of the absorbing ga s. They are consistent with the find- 
ings of IKacprzak et alj d2010al) that the optical depth and velocity 
spread of the gas, relative to the galaxy systemic velocity, increases 
wi th disk inclination. Interesting ly, we do not find a trend suggested 
by ICharlton & Churchilll J 19961) for increasing W rs (normalized by 
the Tully Fisher velocity, computed from L B ) with increasing incli- 
nation (0.2cr). A simple galaxy disk model is not supported by the 
position angle data, which are not necessarily consistent with the 
co-planer geometry scenario. If the absorbing gas were in a disk, 
then the radial distance out on the disk mid-plane probed by the 

line of sight is given by r = D V 1 + tan 2 i sin 2 8. Correlation tests 
between the absorption properties and r might be expected to have 
higher significance (less scatter) than those with D (even allowing 
for all disks not being equal). For our sample, the absorption prop- 
erties exhibit uncorrelated scatter with r (no test had a significance 
above 2a). Thus, at best, we find the data suggestive that most Mgn 
absorbers have a co-planer geometry, but are not n ecessarily disks 
(also s ee the geometric and kinematic models of IKacprzak et al .1 
l2010al) . Since there is a substantial scatter in the absorption proper- 
ties, there is some difficulty in interpreting these trends. Since the 
distribution of inclinations for our sample is consistent with that of 
a random sample of galaxies, our sample has a larger number of 
higher inclination galaxies than lower inclination galaxies [18/40 
galaxies have i > 60° and only 5/22 have i < 30°). With a paucity 
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log(D) (kpc) log(D)-logO) 

Figure 6. — (a) The 3.2cr anti-correlation between W r (2796) and impact parameter D. Note the large scatter in the correlation. The dark line shows a 
maximum-likelihood fit to the data with log( W,-(2796)) = 1.228 - 0.8841og(Z>). The points are color-coded as a function of galaxy inclination angle, i. It is 
apparent that galaxies with high, mid, and low inclinations populate different regions of the plot. For a fixed D, high inclination galaxies tend to have higher 
W, (2796) than low inclination galaxies. This inclination gradient in the W,-(2796) direction suggests that inclination plays a secondary but significant role in 
determining the W,-(2796) of an absorption system. — (b) A correction for inclination is applied to the W,-(2796) and D anti-correlation. The scatter is now 
reduced and results in a 4.3cr correlation between W r (2796) and i/D. Note the different inclinations are now overlapping compared to (a). The dark line shows 
a maximum-likelihood fit to the data with log(W,.(2796)) = -0.334 - 1.3101og(D/i). 



of low inclination (predominantly face on) galaxies, any possible 
scatter in the absorption properties for low inclination may not be 
well represented in our sample. However, there is no clear or evi- 
dent bias in our sample that could explain why the lower inclination 
galaxies have lower strength absorption properties. 



3.7 Normalization by Impact Parameter 

As discussed in Section 13.31 the anti-correlations between 
W r (2796) and D and between W r (2803) and D, and decreasing 
trends with the remaining absorption properties, all corroborate a 
physical picture in which the column densities diminish with in- 
creasing D. This suggests that impact parameter has a strong in- 
fluence on the absorption properties and that calibrating out im- 
pact parameter for the correlation tests between galaxy properties 
and absorption properties may uncover otherwise diluted physi- 
cally motivated trends or correlations. From the point of view of 
non-parametric rank correlation estimators, multiplying the absorp- 
tion properties by D is equivalent to dividing the galaxy properties 
by D. We adopt the latter approach. Since one of the strongest in- 
dicators of a dependence of the absorption on the galaxy properties 
is orientation, we focused on the galaxy inclination^ 

Following normalization of i by D, we find that N vp , N c i, r, 
N a , W, (2803), and W r (2796) are all correlated with the normalized 
inclination, i/D, at a greater than 3.2<x significance. The strongest 
correlation is between 1^, (2796) and i/D, which has a one part in 
1(P probability of being consistent with no correlation (4.3cr). The 



6 We performed the correlation tests on all galaxy properties, but find that, 
unless discussed in the text, there were no other statistically significant (i.e., 
greater than 3cr) correlations when normalizing by D. 

7 We have removed galaxy Q2206-199G2 from the rank correlation tests 
when normalizing by D because its value of inclination and errors, (' = 
2.9+' 7 g 4 and D = 105.2kpc, provides no constraints on the ratio of log(D/i). 
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Figure 7. The galaxy inclination as a function of the equivalent width resid- 
uals computed from the fit between W r (2796) versus D (Figure[6^). We find 
a 2.6<t correlation between the residuals (scatter) and the galaxy inclina- 
tion. This correlation further demonstrates the significance in the reduction 
of scatter between W, (2796) versus D when galaxy inclination is taken into 
account. 



results are presented in Table[6] We note that D does not appear to 
correlate with ( (only at 1.5cr). 

In Figure [6^, we show the 3.2<x anti-correlation between 
W, (2796) and D. As we have discussed in Section [3"3l there is large 
scatter observed in the W,-(2796) and D anti-correlation. In Fig- 
ure^, we also color-code the data points as a function of galaxy in- 
clination. It is clearly noticeable that high, mid, and low galaxy in- 
clinations populate different regions of the plot and are offset from 
each other as a function of W, (2796) (i.e., for a fixed D, high incli- 
nation galaxies tend to have higher W, (2796) than low inclination 
galaxies). The galaxy inclination gradient in the W,(2796) direc- 
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tion suggests that inclination plays a secondary role but significant 
determining the W r (2796) of an absorption system. 

In Figure [6p, we apply a galaxy inclination correction. We 
note that our discussions of rank correlation tests are in terms of the 
galaxy inclination normalized by the impact parameter, i/D, how- 
ever, in Figure|6p we chose to plot D/i in order to aid the reader to 
directly compare the scatter from both panels. Note that the scatter 
is significantly reduced and that there is now no clear inclination 
gradient seen in the distribution of data. A maximum-likelihood fit 
to the iy, (2796) and i/D correlation is provided in Table|7j In Fig- 
ure|7jwe show that if we fit the anti-correlation between W,-(2796) 
versus D (Figure^), then the equivalent width residuals computed 
from the fit correlate with galaxy inclination at 2.6cr, further sup- 
porting the significance in the reduction of scatter in Figure |6p. 

We also apply bootstrap statistics to test the probability that a 
random distribution of inclination values could produce the incli- 
nation gradient as a function of W, (2796) seen in Figure^, which 
results in the tightening of the anti-correlation seen in Figure[6p. We 
used our 40 data points of W, (2796) and D pairs since they exhibit 
the expected correlation shown in Figure^. We then took our mea- 
sured values of i and randomly reassigned them to the W, (2796) 
and D pairs, thereby creating a new bootstrap realization of the 
sample. We then repeated this process 1 x 10 6 times and computed 
the Kendall rank correlation test for each realization of the sample. 
We find that the probability of this correlation occurring by chance 
due to random selection of galaxy inclinations is P = 0.00040. 
Thus, it is unlikely that this correlation is due to random chance at 
the 3.54<r significance level. 

Thus, we find that the distribution of halo gas absorption 
strengths has both an impact parameter and inclination depen- 
dence, which suggests that the gas in the halo is not spherically 
distributed. That is, these results suggest that Mgn halo gas has a 
co-planer geometry and is coupled to the inclination of the galaxy 
disk. The correlation is opposite to that expected for wind models. 

We also examined if multiples of the galaxy bulge and disk 
scale lengths, D/ri, and D/rj, and half-light radii, n,, might further 
reduce the scatter in the above correlations by using the quantities 
i/(D/r d ), etc. However, we found that multiples of the scale lengths 
resulted in greater scatter and no statistically significant correla- 
tions with absorption properties. We also normalized the galaxy 
position angle by D and re-examine the scatter in the absorption 
properties with 6. We found a slight increase in the significance 
levels from 1.4<x and below to ai 3. Oct and below. However, these 
significance levels remain lower than the significance level of the 
anti-correlation between W r (2796) and D. As such, the increased 
significance levels are primarily driven by the normalization of D. 
The upshot is that, if the selection of galaxies by Mgn absorption 
yields a preferred position angle-impact parameter range in rela- 
tion to the absorption strength properties, it has a much larger scat- 
ter than does the preferred inclination-impact parameter range for 
a given range of absorption properties. 



4 DISCUSSION 

Several studies have shown that the morphologies of intermedi- 
ate redshift (0.3 < z < 1.0) Mgn absorption-selected galaxies 
appear to be qualitativ e ly similar to those of "typical" local field 
galaxies JSteidelll 19981 ; IChen et alj|20ol IChen & Lanzettall2003l ; 
IKacprzak et alj|2007l) . Using GIM2D to quantify the morphologi- 
cal parameters of 0.3 < z < 1.0 Mgn absorption-selected galaxies, 
we find that their bulge and disk size distributions are also similar to 



those found for field galaxies at similar redshift. Our results quanti- 
tatively show that galaxies selected by Mgn absorption are, broadly 
speaking, "typical" galaxies. Nevertheless IKacprzak et al.l J2007h 
has shown that Mgn absorption-selected galaxies appear to have 
higher morphological asymmetries than field galaxies. If the latter 
is true, then accurately modeling absorption-selected galaxies with 
smooth light profiles may be more difficult given their higher level 
of morphological asymmetries. The difficulties may induce more 
scatter in the distribution of morphological parameters thereby po- 
tentially mitigating the significance levels of any correlations be- 
tween galaxy properties and absorption properties. 

It remains strongly debated whether Mgn absorption systems 
arise from star formation driven winds or an array of structures 
such as tidal streams, satellites, filaments, etc. Here we argue that 
our sample of galaxies, 75% of which have W,-(2796) < 1.0 A, 
does not support a scenario in which winds are the predominant 
mechanism producing the absorption. Overall, it is likely that an 
admixture of these processes contribute to the observed Mgn ab- 
sorption systems; below, we will argue that winds may dominate 
the high equivalent width regime, whereas other processes domi- 
nate the lower equivalent width regime. 

For reference, in Table |7J we list the known quantified corre- 
lations between the Mgn absorption properties and the host galaxy 
properties, which we will be discussing in the remainder of this 
section. The W r (2196) range used in each study is also listed. 

Studies of strong Mgn systems (i.e., W,-(2796) > 1 A), 
where host galaxies typically have D < 25 kpc, have shown 
that they are likely produced by winds from h i gh sta r-forming 
galaxies (e.g.. |Prochter et al ] |2006h . lBouche et alj Hoolih report an 
anti-correlation between W, (2796) and the amplitude of the cross- 
correlation betw een luminous red galaxies (LRGs) and M gn ab- 
sorbers (also see iGauthier et al]|2009l : iLundgren et"ai1l2009h . They 
conclude that stronger absorption systems are not produced by viri- 
alized gas within g alaxy halos but that they originate in supernovae- 
dr iven winds. T he Bouche et_alj results are consistent with those 
of IZibetti et alj d2007» . who find that stronger absorption systems 
with W r ( 2796) > 1.2 A are associated with blue star-forming 
galaxies. Mena rd et alj J2009h showed that for stronger systems, 
W,.(2796) correlates with the associat ed [On] luminosity, an e sti- 
mator of star formation rates (also see iNoterdaeme et alKOld) . In 
a detailed study of t wo ultra- strong Mgn absorption systems, with 
W r (2796) > 3.6 A. iNestor et all feoid) reported nearby galaxies 
with star formatio n rates exceeding ~ 10-90 M Q yr~' (also see 
INestor et al. Eool) . These results support the idea that winds pro- 
duced by star-forming regions and/or supernovae are responsible 
for a large proportion of the strong absorption systems (> 1 A) 

detected in the halos of galaxies. 

Our study, and oth ers (e.g., IChurchill. Steidel. & Vogtl 

! 19961: IChurchill & Vogtl 1200 ll ; ISteidel et alj 20021 : 
Churchill. Vogt. & Charltor] l200ll; IChen & Tinked 120081 : 
Chen et alj l2010allbl : IKacprzak et al] l2010al) . use higher reso- 



lution spectra that are sensitive to lower equivalents widths. The 
general conclusions from these moderate-to-low equivalent width 
studies all point to other sources for producing the observed Mgn 
absorption. 

Using a simila r sample to the one presented here, 
IKacprzak etail d2007t) found a correlation between galaxy asym- 
metry, normalized by D, and W, (2796) suggesting that galaxy mi- 
nor mergers and harassments may be producing the absorption de- 
tected in halos. Their correlation strengthens when W,(2796) > 
1.4 A systems were removed, suggesting that interactions (or pro- 
cesses that gently perturb galaxy morphology) may dominate the 
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Table 7. Currently known quantified correlations between the Mgn absorption properties and the host galaxy properties: (1) The correlation between W r (2796) 
and galaxy inclination normalized by impact parameter, i/D; (2) The anti-correlation between D and Mg with W,(2796); (3) The correlation between [On] 
luminosity surface density (ZL^Q n j) and W r (2796); (4) The correlat ion between galaxy asymmetry (A), normalized by D, and W,-(2796); (5) The anti- 
correlation between the host galaxy halo mass (M; m ; D ) and VW(2796). lBouche et al.1 |2006I) did not provided parametrized relationship for their correlation, 
however we used the binned data from their Table 3 and applied a least least squares fit. 



Properties 


Redshift 
Range 


W,(2796) 
Range (A) 


(Anti-) 
Correlation 


Signif- 
icance 


Parametrized Relationship 


Reference 


i/D & W r (2796) 


0.3-1.0 


0.03-2.9 


Correlation 


4.3tr 


log(W r (2796)) = -0.334 - 1.3101og(D/0 


This paper. 


D + M B & W,-(2796) 


0.1-0.5 


0.1-2.4 


Anti-correlation 


7<r 


log(W r (2796)) = Alog(Z)) - B(M B - M* B ) + C 
A = -1.93 ±0.11, B = -0.27 ±0.02 
C = 2.51 ±0.16 


Chen et al. (2010a) 


ZL [Qll] & H/, (2796) 


0.4-1.3 


0.7-6.0 


Correlation 


15o- 


SL [0n] =A(W,(2796)/1A) ' 

A = 1.48 ± 0.18 X 10 37 erg s _1 kpc 2 

a = 1.75 ±0.11 


Menard et al. (2009) 


A/D & W r (2796) 


0.3-1.0 


0.03-2.9 


Correlation 


3.3tr 


A/D = 2.48 x 10- 3 W,-(2796) + 8.15 x 10~ 4 


Kacprzak et al. (2007) 


M hah & W,-(2796) 


0.4-0.8 


0.5-5.0 


Anti-correlation 




\og(M halo ) = 13-0.6W r (2796) 


Bouche et al. (2006) 



Mgn absorption profiles in the lower equivalent width regime. For 
our sample, we find trends that weaker absorption properties are 
associated with redder galaxies; but there is no suggested trend di- 
rectly wit h ty r (2796) a s found for higher W, (2796) systems as re- 
ported by IZibetti et alj d2007l) . The lack of statistical significance 
may be due to the small number of galaxies in our sample. For a 
sample of 71 absorption sel ected galaxies with a W r (2796) distri- 
bution similar to our sample. lChen et alJoOlOal) also did not find a 
significant correlat ion between !V r (2796) and gala xy color. The dis- 
crepancy between Izibetti et"al] and our work and lChen et alj may 
indicate that there is a fundamental difference in the galaxies se- 
lected by stronger and weaker Mgn absorption systems; weaker 
systems are possibly not correlated with host galaxy star-formation 
rates or the luminosity-weighted stacking procedure is not fully un- 
derstood. However, this correlation may have a subst antial scatter, 
so that large samples, on the order of those used by IZibetti et all 
are required to obtain statistical significance. 

The Mgn equivalent widths, and the other measures of the 
Mgn column density, decrease with D, indicating that halos exhibit 
a natural decrease in gas column density with increasing projected 
galactocentric distance. This anti-correlation alone cannot be lever- 
aged to differentiate between wind scenarios and other Mgn gas 
producing mechanisms; however, we would expect a difference in 
the halo gas geometry for the competing scenarios (see § 13.61 and 
§ |3,7| l. Our sample exhibits trends for increasing absorption prop- 
erties (and scatter) with higher inclination galaxies. Accounting for 
the decreasing gas column density (and its scatter) with increasing 
D reduces the scatter and yields a strong correlation (4.3cr). At the 
very least, this suggests that the combined effect of inclination and 
impact parameter is such that weaker absorption is found at larger 
D and smaller i (far from face-on galaxies). As i increases for a 
given D, we see W,(2796) increases, or as D decreases for a given 
i, we see W, (2796) increases. Given that the velocity widths, W vs , 
do not follow this behavior, the data favor an increasing column 
density (path length) for increasing inclination at a fixed D. Such 
behavior is expected for co-planer geometry with kinematics that 
are not strictly coupled to disk rotation, but is not what is expected 
for a wind geometry. For a simple wind scenario, an anti-correlation 
between W, (2796) and i/D would be expected. These trends are not 
evident in our data. 

Recent cosmological SPH simulations o flStewart et alj J201 ll) 



have shown that gas-rich mergers and cold-flow streams can pro- 
duce a circum-galactic cool gas component that predominately in- 
falls towards the host galaxy disk. They note that this gas accretion 
should be observed as a relative host-galaxy/halo-gas velocity off- 
sets of ~ 100 km s" 1 . These offsets have already bee n directly ob- 
served at z ~ 0- 1 - 1 -0 using Mgn a bsorption systems dSteidel et al.1 
|2002| ; iKacprzak et ai]|2010j, l201ll) . In such a scenario one would 
expect a correlation between host galaxy inclination and absorption 
strength, column density, etc., if the circum-galactic component 
was a pure extension of the hos t galaxy disk. According to these 
simulations. IStewart et alj d201 ll) found that in most cases the ac- 
creting gas co-rotates with the central disk in the form of a warped 
extended cold flow disk, and the observed velocity offset are in the 
same direction as galaxy rotation. These models support the corre- 
lations between the absorption properties and i/D found here, and 
the warps observed in the simulations may further explain the scat- 
ter seen in these correlations. 

The relative gas-galaxy kinema tics also differ betw een low 
and high equivalent width systems. iBond et af I l l2001bl) studied 
the kinematics W r (2796) > 1.8 A systems and found character- 
istics of wind driven gas (double quasi-symmetric absorption over 
~ 400 km s _1 , not characteristic of Mgn in DLAs). On the other 
hand, systems with W,(2796) < 1 A are characterized by velocity 
widths no greater than ~ 50 km s" 1 accompanied by one to a few 
very weak a nd narrow components a n d are not suggestive of wind 
kinematics dChurchill & Vogtl l200ll ; IChurchill. Vogt. & Charltor] 
l200ll). For predominantly W r (2796) < 1 A systems, IChen et al. 
( 120 lOaf) found that the velocity differences between the Mgn ab- 
sorption and the host galaxies are roughly 16 km s" 1 with a dis- 
persion of 137 km s~'. These velocity offsets are much lower than 
expected for wind driven outflows. 

In fact, a direct kinematic comparison of six galaxies and their 
M gn absorptio n kinematics, with velocity offsets similar to those 
of IChen et all were shown to be fairly well described as having 
lagging disk-like kinematics, th ough not all of the velocity spread 
could be successfully modeled dSteidel et alj|2002l) . However, the 
discrepancies were on the or der of 30 km s~ ' in the most extreme 
cases. For a similar sample, IKacprzak et al.1 d2010al) also showed 
that, even if the Mgn absorption tends to reside fully to one side 
of the galaxy systemic velocity and aligns with one arm of the ro- 
tation curve, not all the absorbing gas kinematics can be explained 
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by a co-rotating halo model. Using cosmological simulations, they 
further showed that even if the majority of the simulated Mgn ab- 
sorption arises is an array of structures, such as filaments and tidal 
streams, the halo gas often has velocities consistent with the galaxy 
rotation velocity. The line of sight gas motions along these struc- 
tures in the simulatio ns reproduc e a velo city distribution consistent 
with that reported bv lChen et ail l l2010ij) . 

Finally, it is also interesting to note the equivalent width red- 
shift path density evolution of Mgn absorbers: higher equivalent 
width systems evolve in that there are fewer high W,(2796) sys- 
tems per unit redshift at low redshifts than at higher redshifts. 
The evolution is more pronounced for the highest W, (2796) sys- 
tems, i.e., those with 2 A and above. As lower W r (2796) sys- 
tems are examined, the evolution weakens such that for all systems 
with W, (2796) > 0.3 A the absorber population is consistent with 
the no-evolution expecta t ion for the presentl y accepted cosmology 
dSteidel & Sargentll 19921 : iNestor et alj|2005h . Assuming winds are 
a transient phenomenon in most galaxies, and given that the global 
star formation rate of the universe decreases toward low redshifts 
dMadau et al.ll 1998b . it would seem reasonable that, globally, the 
higher equivalent width systems have some causal connection to 
star formation rates. On the other hand, the fact that lower equiva- 
lent width systems do not evolve with redshift would indicate that 
the processes producing lower W, (2796), though varied, are fairly 
constant with cosmic time. Though major merger rate s are expected 
to de crease strongly, oc (1 + z) 2A , with redshift (cf.. IStewart et al.l 
1 20091) . it is not clear that other scenarios such as minor mergers 
and/or accretion of the IGM (filaments, etc.) evolve as rapidly. It 
could be that the lack of redshift evolution of lower W, (2796) ab- 
sorbers is linked to the latter, possibly more ubiquitous processes. 

The different physical characteristics and evolution of the pop- 
ulations of high and low W,-(2796) Mgn systems suggest different 
physical mechanisms giving rise to each population. The correla- 
tion between W,(2796) and i/D for our lower W,(2796) sample 
indicates that winds do not dominate the lower equivalent width 
regime. It appears that galaxy inclination plays a strong role in de- 
termining the optical depth of the halo gas once the decreasing col- 
umn density with increasing impact parameters is taken into ac- 
count. We find that edge-on systems are likely to produce higher 
optical depth absorption systems, and lower optical depth absorp- 
tion systems are produced by face-on galaxies: this is the opposite 
effect expected for a wind scenario. These results support a picture 
where the Mgn absorption arises in structures that are relatively 
co-planer to the host galaxy disk. Such structures might include ac- 
creting filaments or tidal streams from minor mergers in the galaxy 
plane, and disk warps. T hese results are consi stent with recent cos- 
mological simulations o flStewart et alj d201 ll) who find that the ac- 
cretion of cool gas via filaments and gas rich merges does in fact 
form a stable disk that supplies gas and angular momentum to the 
host galaxy. 



5 CONCLUSIONS 

We have performed a detailed study of a sample of 40 Mgn 
absorption-selected galaxies between 0.3 < z < 1.0. The galax- 
ies have B-band absolute magnitudes that range between -18.6 < 



measurements from the literature. The galaxies are at projected sep- 
arations of 12 < D < 115 kpc from the quasar line-of-sight. 

We have used GIM2D to model WFPC-2/ffiT images and ex- 
tract quantified morphological parameters for 40 Mgn absorption- 
selected galaxies. These parameters include: the bulge-to-total frac- 
tion (B/T), the bulge semi-major axis effective radius (r b ), the 
bulge ellipticity (e/,), the bulge position angle (Qt), the bulge Sersic 
index (n), the semi-major axis disk scale length (r rf ), the disk incli- 
nation (0, and the disk position angle (0 d ). These properties help 
us further compare and quantify the nature of absorption-selected 
galaxies. 

Furthermore, we have extracted absorption parameters from 
the quasar spectra as well as from Voigt profiles fits to the absorp- 
tion systems. We have measured the optical depth weighted mean 
Mgn /12796 absorption redshift, the rest-frame equivalent width 
W r (2796), the doublet ratio, the number of clouds, the Voigt pro- 
file fitted system column density, the Mgn optical depth, the AOD 
derived column density, the mean velocity, the velocity spread, and 
the velocity asymmetry. 

In order to explore possible connections between the Mgn ab- 
sorption properties and the galaxy morphological properties, we 
have performed non-parametric Spearman and Kendall rank cor- 
relation tests. Our mains results can be summarized as follows: 



(i) Mgn host galaxies appear to be similar to those at low red- 
shift and have a wide range of morphologies and colors. With 27 
late-type galaxies and 13 early-type galaxies, their populations are 
consistent with the distribution of galaxy morphologies found in the 
field environment. They have a mean disk scale length of 3.8 kpc 
and a mean bulge scale length 2.5 kpc, which are comparable those 
of the Milky Way. The disk scale lengths and bulge effective radii 
distributions of the sample are similar to those of field galaxies 
obtained from the DEEP survey, at similar redshifts. The mean 
W / ,.(2796) for the early-type galaxies is 1.0 A and slightly lower for 
late-type galaxies at 0.7 A, where both types have a large spread 
in equivalent width, although late-type galaxies dominate for ab- 
sorption systems with W r (2796) < 0.3 A and tend to also be at 
D > 40 kpc. 

(ii) We find a 3.2<x anti-correlation between D and Mgn equiv- 
alent width. There is large a scatter in the distribution, suggesting 
that D is not the only physical parameter affecting the distribution 
and quantity of halo gas. There are no other absorption properties 
that scale with D above 3<r. However, all our other measures of the 
gas absorption strength show decreasing trends with D. Taken to- 
gether, the anti-correlation trends show a consistent signature that, 
on average, Mgn absorbing gas becomes optically thinner with in- 
creasing projected galactocentric distance from the host galaxy and 
is likely the dominant behavior of the halo gas. 

(iii) We find only weakly significant trends (2 - 2.2cr) between 
color and the absorption properties, suggesting that larger absorp- 
tion quantities tend to be associated with the bluer gal axies. We do 
not rep roduce the B - K and 1^(2796 ) correlation of Zibetti et"al] 
d2007t) . Our results are consistent with lChen et al.1 ( feOlOal) and we 



-23.8 and are associated with Mgn absorption systems 



with rest-frame equivalent widths that range between 0.03 < 
W,.(2796) < 2.9 A. The Mgn absorption profiles were obtained 
from HIRES/Keck and UVES/VLT quasar spectra; we did not have 
high resolution spectra for six systems and use equivalent width 



conc lude that our samples probe a lower equivalent width range 
then lZibetti et all and therefore are likely probing different mech- 
anisms producing Mgn absorption. We find less than a l<x connec- 
tion between M B and the absorption properties, implying that for 
our luminosity and equivalent width range the Mgn absorption is 
not strongly dependent on galaxy luminosity. 

(iv) We find no correlation between 6 and the absorption prop- 
erties (the highest significance is 1.4cr); the absorption properties 
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are consistent with a random distribution as a function of position 
angle. 

(v) By accounting for the decreasing gas column density (and its 
scatter) with increasing D, the correlation with i/D and W,-(2796) 
increases to 4.3<r significance level. Also, following normalization 
of ( by D, we find that N a , r, N ch N vp and VK,(2803) are all corre- 
lated with i/D greater than 3.2tr level of significance. Overall, these 
results suggest that the Mgn gas has co-planer geometry, but is not 
necessarily disk-like, that is coupled to the galaxy inclination. This 
results do not support Mgn absorption produce by sta r-burst driven 
winds . These results are consistent with the models o flStewart et al.l 

feouh . 

(vi) We do not find any other correlations above 2. Oct between 
the remaining galaxy properties, such as r^, n,, r e , B/T, n and the 
absorption properties. Thus, the galaxy optical size and shape does 
not appear to be an important factor in determining the amount of 
gas within a galaxy halo. 

We find several interesting connections between the galaxy 
morphological properties and the Mgn absorbing gas. Although re- 
cent evidence suggests that high equivalent width systems are pro- 
duced predominantly by winds in star-forming galaxies, the results 
of this paper do not support such an explanation for weaker sys- 
tems (W,(2796) < 1 A). The correlation between the inclination 
of the galaxy disk and the halo gas absorption strength is suggests 
that the the halo gas resides in a co-planer distribution. It is plausi- 
ble that the absorbing gas arises from tidal streams, satellites, fila- 
ments, etc. which tend to have more-or-less co-planer distributions. 
This correlation could not be explained in the wind scenario. 

One of the few ways to differentiate between winds and other 
sources replenishing the halo gas is by studying multi-phase gas 
sensitive to different density and temperatures. These questions ex- 
press the necessity of UV spectrographs like COS and STIS. 
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APPENDIX A: APPENDIX A: SUPPORTING 
INFORMATION 

Figure Al is the complete version of FigureQ] It will appear in the 
online version of this paper, not in the printed version. 
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Figure Al. — (far-left) The HIRES/Keck or UVES/VLT quasar spectra of the Mgn /12796 absolution feature are shown alongside the associated absorbing 
galaxy on the right. The Mgn /12796 optical depth weight mean absorption redshift is the zeropoint of the velocity scale. The tick marks indicate the number 
of Voigt profile components and the red curve indicates the fit to the data. We do not have HIRES or UVES data available for six galaxies, (left) WFPC-2/HST 
images of galaxies selected by Mgn absorption. The images are 10 times larger than the 1.5cr isophotal area. — (center) The GIM2D models of the galaxies, 
which provide quantified morphological parameters. A scale of one arcsecond is indicated on each image along with the physical scale computed at the Mgn 
absorption redshift. — (right) The residual images from the models, showing quality of the fit and the underlying structure and morphological perturbations of 
the galaxies. The encircled arrow provides the direction to the quasar (galaxy-quasar orientation). The cardinal directions are also shown and the quasar name 
and redshift of Mgn absorption is provided under each set of galaxy WFPC-2, model and residual image. 
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Figure Al. — continued 
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